Exactly solvable analogy of small- world networks 
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We present an exact description of a crossover between two different regimes of simple analogies 
of small-world networks. Each of the sites chosen with a probability p from n sites of an ordered 
system defined on a circle is connected to all other sites selected in such a way. Every link is of a 
unit length. Thus, while p changes from to 1, an averaged shortest distance between a pair of sites 
changes from I ~ n to I = 1. We find the distribution of the shortest distances P(£) and obtain 
a scaling form of £(p,n). In spite of the simplicity of the models under consideration, the results 
appear to be surprisingly close to those obtained numerically for usual small- world networks. 
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A sharp interest on small-world networks began last 
year, when a beautiful paper of Watts and Strogatz 
was published. A number of publications on the topic is 
rapidly increasing | 0. One of versions of the model 
proposed by Watts and Strogatz looks as the following 
(see ref. ||). Let n sites be arranged on a circle. Each 
of several neighbours are connected by unit length links. 
Then each of those links is rewired with a probability p 
to one of other randomly chosen sites. Thus, one stud- 
ies a crossover from an ordered structure at p = to a 
random graph [fLBj. 

The main quantity which is considered for these mod- 
els is the shortest distance between two sites of a network 
(that is a quantity which is much more difficult to cal- 
culate than connectivity, for instance). When p = 0, the 
averaged over all pairs of sites shortest distance £ is equal 
to n/4, and l~ logn if p = 1 Q. 

The central result obtained by Watts and Strogatz 
is very simple: even addition of a single shortcut to a 
large system decreases I noticeably, so that £ is a quan- 
tity extremely sensitive to such "global" changes in a 
network as addition of a shortcut between far sites. 

The main problem of the small network theory is to de- 
scribe the crossover between the p = and p — 1 regimes 



of a finite size network and to find a scaling form 19 for 
£(p,n). One should emphasize that in spite of a lot of 
effort to understand properties of small-world networks, 
the only exactly proved statement for £ was the following: 
the scaling variable of the crossover under consideration 
is pn - e.g., see ||. The reason of the failure is the dif- 
ficulty of calculation of I. In the present communication 
we consider very simple models with such a crossover be- 
tween states with two different values of £ for which those 
calculations may be made exactly. We shall present the 
scaling form of £ and we shall describe P(£) analytically. 
The answers are surprisingly close to the corresponding 
results obtained numerically for previously studied much 
more complex small- world networks. 

Let us introduce the models. Let n sites be arranged 



on a circle. Each of them is connected with two near- 
est neigbhours by links of unit length. We shall start 
our consideration from the most simple for consideration 
model with all the bonds on the circle oriented in one 
direction as shown in Fig. 1(a). Afterwards we shall 
study the model with all undirected links [see Fig. 1(b)]. 
We add a new central point to a network and connect it 
with any other site chosen with the probability p by an 
undirected link of length equal to 1/2. In fact, it means 
that we select some sites with the probability p and then 
connect all of them by undirected unit length links. 

Let us find the distribution P(£) of the shortest paths 
between pairs of sites. The calculations are most simple 
for the model with directed bonds showed in Fig. 1(a). 
In the present short communication, we use, maybe, the 
simplest way to obtain the distribution. At first, we 
shall obtain the probability P(£,k), i.e. the probabil- 
ity that the shortest path from one site to another is £ 
when the distance between them counted along the ring 

is k: J2e=i p (^ k ) = !• We sha11 nnd p ( e < k > k ) for 
several small values of k and shall demonstrate that it 
can easily be described by a general formula for all £ < k 
and k. 

To find P(£,k) for the model shown in Fig. 1(a), one 
has to consider all possible configurations of links con- 
necting the center and sites i = 0, 1, 2, . . . , k and to cal- 
culate statistical weights for the shortest path to be equal 
to a particular value £. In such a way, one may easily ob- 
tain the following expressions for small £ and k: 
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The next expressions for P(£, k) show the same trend. 
Thus, the general form of P(£, k) is 
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The distribution of shortest distances between pairs of 
sites for the model under consideration looks as the fol- 
lowing: 

n — 1 -, n — 1 

p w - — E ^. fc ) = — E p (^> fc )- (3) 

fe=i k=e 
Inserting Eqs. (2) into Eq. (3) we get immediately 

P(£) = — !— - [l + (* - l)p + *(n - 1 - l)p 2 ] (1 - p)* -1 . 
n — 1 

(4) 

If p 0, P(£) -> l/(n - 1) for any £ < n - 1. If p ->■ 1, 
P(l) — > 5« i , where is the Kronecker symbol. 
The average shortest distance, £ may be written as 



^£P(£), 



(5) 



so using the previous distribution one gets 
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Thus, £[p — > 0) — > n/2 and £(p — > 1) — > 1. The maximum 
value of the distribution, P(£) is at the following point, 
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Again, £ ma x{p -> 0) -> n/2 and l max {p -> 1) — > 1. 

The Eqs. (4) and (6) are our main results for the 
defined model. To obtain a scaling description of the 
crossover region, one has to pass to the limits n — > oo 
and p — > while the quantities p = pn and z = £/n are 
fixed. In this limit, from Eq. (4) one obtains a continuous 
distribution Q(z,p) of z: 



nP{£) = Q(z, p) = [1 + pz + p 2 z{\ - z)\ e- pz 



(8) 



< z < 1. In Fig. (2) we present Q(z,p) for several val- 
ues of p as a function of shortest distance z. The average 
value of the shortest distance may be found from Eq. (6) 
or (8): 
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The asymptotic behavior of z is 1/2 — p 2 /24 (p — ► 0) 
and 2/p (p ^> 1) and of z max is 1/3 — 2p/27 (p — > 0) 
and 1/p (p> 1). From Eqs. (||) and (||), one sees that 

(z - z) 2 ~ 2/p 2 when p > 1, so J {z - zfjz — * l/\/2 
in that limit. 

In Figure 3, we plot ^ and £ ma x vs p. Note, that the 
plot of £/n resembles strikingly the results obtained nu- 
merically (e.g., compare with [|l]J^||) and by series expan- 
sions j fL2[ for ordinary small- world networks, although we 
consider much more simple network. The analytical re- 
sult for the distribution Q{z) presented in Fig. 2 looks 
very similar to the corresponding numerical one obtained 
by Barrat and Weight ||] for a standard small-world net- 
work. 

The model we have studied above looks in some sense 
rather artificial because of the directed links on the circle. 
Now we shall consider a more naturally looking model 
without directed links [see Fig. 1(b)]. We can proceed in 
the same way as before to obtain the following probabili- 
ties for small k (although, the procedure is more tedious 
in this case, since now to obtain P{£, k) one has to count 
over all configurations of links connecting the center and 
sites -k + 2, -k + 3, . . . , 2k - 3, 2k - 2), 
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(One may check that Eqs. ( |Tl| ) hold for any possible n.) 
As before, these relations may be written in the following 
general form, 

P(l,l) = l, 
P(£=l,k)=p 2 , 
P(2<1< k,k) = P 2 (l-p) 2l - 4 (2- P )(2l-2-£p), 
P(£ = k,k) = (l- P ) 2k ~ 4 [l + (2k - A)p - (k - 1) P 2 } 
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For an odd 71, for instance, one can write 

„ (n-l)/2 



fe=l 



Substituting Eq. @ into Eq. (||) one gets 



P(l) 
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and 



P(i > 2) = [2 + 4(£ - 2)p + 2(£ - l)(2n -At- i)p 2 - 
2(21 - l)(n -2*- l)p 3 + £(n - 2* - l)p 4 P ^ . 

71—1 

This expression for P(£) looks slightly different from 
Eq. (|]) for the previous model. But when we take 
the scaling limits, we find the continuous distribution 
nP(£) = Q'(z, p) — 2Q(2z, p), which differs from the cor- 
responding distribution Q(z,p) for the model with ori- 
ented links [see Eq. (0)] only by a scaling factor of two: 
z — ► 2z (now < z < 1/2). Therefore, after that trivial 
change of a scale, Eqs. (||) and (|i"tj| ) and plots in Fig. 3 
are also valid in the present case. In fact, the results for 
both models under consideration turns to be the same. 

In summary, we have studied two networks which are 
simple analogies of small- world networks. We have ob- 
tained explicit expressions for the distribution P(£) of 
the shortest paths between pairs of sites and for the av- 
erage shortest path I, The scaling functions look similar 
to the corresponding ones obtained numerically for more 
complex small- world networks [|]||pj. The studied net- 
works (especially, the second one) model a rather real 
for our world situation in which all far connections occur 
through some common center. 

Note that we have used here the simplest but not 
the most rigorous way to find an exact solution of the 
models under consideration. In fact, we have demon- 
strated above how to guess the solution, so the problem of 
searching for explicit solutions for small-world networks 
remains open. 



One should notice finally, that the above considered 
models differ from the original small-world networks [jl]j|] 
in two aspects. One of them is evident: the studied mod- 
els are of mean-field nature since all the selected sites are 
connected together by shortcuts. Another aspect is that, 
although, the original small- world networks are in an or- 
dered state for p = and they are random graphs for 
p = 1, in our case, the networks are ordered for both 
p = 0, 1. Nevertheless, as we have shown, the results for 
our simple models resemble strikingly the corresponding 
answers for the original small- world networks. What is 
the reason for such a similarity? First of all, we have 
to say that the obtained crossover behavior seems to be 
independent of the mean-field nature of our models. Sec- 
ond, the scaling regime is observed in the limit p — > 
when all such networks are in the disordered state. Nev- 
ertheless, we have solved analytically the model in which 
the links connected with the central site in Fig. 1 are 
directed chaotically fl20|| , so the crossover is between an 
ordered and disorded states. The obtained results for I 
again look similar to Fig. 3. 
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FIG. 1. Models under consideration, n sites on the circle 
are connected by unit length links. A central site is connected 
to a particular site chosen with the probability p by a link 
of half length. The links connected with the center are not 
directed, a) All links on the circle have a defined orientation, 
b) There are no directed links in the network. 




FIG. 2. The distribution Q(z,p) of normalized shortest 
paths z = l/n between pairs of sites for different values of 
p = pn for the model shown in Fig.l (a) [see Eq. (^)]. Curves 
labeled by numbers from 1 to 6 correspond to p — 0, 2, 5, 8, 
11 and 14. For the model without directed links, the corre- 
sponding distribution equals 2Q(2z,p). 
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FIG. 3. The scaling functions z = I/n and z max = £ m ax/n 
vs. p = pn for the model with directed links. The upper curve 
corresponds to z. For the model without directed links, one 
should change z — > 2z and z max — > 2z max . 



